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Abstract 

Higgs bosons can be searched for in the channels pp/pp — ► ttH + X at the Tevatron and the LHC. We 
have calculated the QCD corrections to these processes in the Standard Model at next-to-leading order. 
The higher-order corrections reduce the renormalization and factorization scale dependence considerably 
and stabilize the theoretical predictions for the cross sections. At the central scale /i = (2mt + Mjj)/2 
the properly defined K factors are slightly below unity for the Tevatron (K ~ 0.8) and slightly above 
unity for the LHC (K - 1.2). 
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Higgs bosons can be searched for in the channels pp/pp — > ttH + X at the Tevatron and the LHC. We have 
calculated the QCD corrections to these processes in the Standard Model at next-to-leading order. The higher- 
order corrections reduce the renormalization and factorization scale dependence considerably and stabilize the 
theoretical predictions for the cross sections. At the central scale /j, — (2m t + Mh)/2 the properly defined K 
factors are slightly below unity for the Tevatron (K ~ 0.8) and slightly above unity for the LHC (K ~ 1.2). 



1. The search for Higgs bosons Q| is one of the most im- 
portant experimental programs in high-energy physics. 
If successful, a crucial step in revealing the mechanism 
for electroweak symmetry breaking and the generation 
of masses for the fundamental particles in the Standard 
Model (SM), electroweak gauge bosons, leptons, and 
quarks, will have been taken. In the near future, the 
search for Higgs bosons will be carried out at hadron col- 
liders, the proton-antiproton collider Tevatron [0 with 
a center-of-mass (CM) energy of 2 TeV, followed by the 
proton-proton collider LHC || with 14 TeV. Analy- 
ses of precision electroweak data H set the focus on 
Mh < 200 GeV as the preferential Higgs mass range in 
the SM, although a firm prediction without escape roads 
is not possible H. 

Various channels can be exploited at hadron colliders 
to search for a Higgs boson in the intermediate mass 
range. Among these channels Higgs radiation off top 
quarks H plays an important role: 



pp/pp — * ttH + X via 



qq, 99 



ttH. 



(1) 



Although the expected rate is low at the Tevatron, 
a sample of a few but very clean events could be ob- 
served for Higgs masses below 140 GeV, while this chan- 
nel becomes very demanding above 0. At the LHC as- 
sociated production of the Higgs boson with top quarks 
is an important search channel for Higgs masses below 
~ 125 GeV. Moreover, analyzing the ttH production rate 
at the LHC can provide information on the top-Higgs 
Yukawa coupling, assuming standard decay branching 
ratios B, before model independent precision measure- 
ments of this coupling will be performed at e + e~ colliders 

!■ 

2. Predictions for the cross sections (M, which are based 
on the leading order (LO), are plagued by considerable 
uncertainties due to the strong dependence on the renor- 
malization and factorization scales, introduced by the 
QCD coupling and the parton densities (see the figures 
below). While estimates of radiative corrections have 
been presented before in the "effective Higgs approxi- 
mation" (EHA) [llO), in this letter we present the first 



complete calculation of the QCD corrections at next-to- 
leading order (NLO), which reduce the spurious scale de- 
pendence significantly and lead to stable predictions for 
the cross sections. 

The Born diagrams, generic examples of which are dis- 
played in Fig. 0(a), are supplemented in NLO by virtual 
gluon-exchange diagrams, Fig. 0(b), running in complex- 
ity up to pentagons; by gluon radiation, Fig. |](c); and 
by parton splitting, Fig. [j](d). The latter two add inco- 
herently to the virtual corrections. 
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FIG. 1. A generic set of diagrams (a) for the Born level, 
(b) for virtual gluon exchange, (c) gluon radiation and (d) 
parton splitting in the subprocesses qq, gg — > ttH etc. 



Dimensional regularization has been adopted for iso- 
lating the ultraviolet, infrared and collinear singularities. 
Renormalization and factorization are performed in the 
MS scheme with the top mass defined on-shell. The top 
quark is decoupled from the running of the strong cou- 
pling a s (/j,). For the evaluation of the pp/pp cross sec- 
tions we have adopted the CTEQ4L and CTEQ4M Q 
parton densities at LO and NLO, corresponding to the 
QCD parameters A^ = 181 MeV and Ap = 202 MeV 
at the one- and two-loop level of a s (/z), respectively. The 
strength of the SM Yukawa coupling is fixed by g U H = 
mt/v, where v = 246 GeV is the vacuum-expectation 
value of the Higgs field, and the top-quark mass is set 
torn* = 174 GeV. 

The most complicated one-loop diagrams are the pen- 
tagons, both analytically and numerically. To calculate 
a five-point integral E^ in D dimensions, the singular- 
ity structure E„-l in D dimensions is determined first. 

" sing 



The singular part E^[ 



(D) 



sing 



is given entirely by three-point 

D) 
sing 



subintegrals. The difference E^ D ' — E^L is finite and 
regularization-scheme independent. Therefore it can be 
calculated in the convenient mass regularization scheme 
in 4 dimensions. The original integral E^ D ' then reads 



E {D) = E (D) 
sing 



E {x 



■ ;D=4) 



7^(mass;D— 4) 
sing 



(2) 



in the limits D -> 4 and mass -> 0. Since £t mass;D=4 ' 
can be expressed in terms of four-point functions |p.2| , 
the D-dimensional five-point integral E*- > is finally re- 
duced to three- and four-point functions. These integrals 
and their tensor structures can be treated according to 
standard methods [fL3|. Numerical instabilities caused 
by vanishing Gram determinants near the phase-space 
boundary, can be controlled by careful extrapolation out 
of the safe inner phase-space domains. [Technical details 
will be presented in a subsequent publication.] 

To extract the singularities of the real part of the 
NLO corrections <j r , a generalization of the dipole sub- 
traction formalism J14| to massive quarks |15| has been 
adopted (see also Ref. |]16| ) . The singularities of the cross 
section cr roal are mapped onto a suitably chosen auxiliary 
cross section a suh which is still simple enough so that the 
singular regions in phase space can be integrated out ana- 



lytically, while the difference a 



real 



a can safely be in- 



tegrated numerically in 4 dimensions. The auxiliary cross 
section c sub can be decomposed into a part crf ub that, de- 
fined on configurations with LO kinematics, cancels the 
soft and collinear singularities of the virtual corrections; 
and a second part u & ^° that includes the singularities 
from initial-state parton splitting, which are absorbed in 
the renormalization of the parton densities. Thus the 
total NLO correction Acr NLO may be written as the sum 



Aa NLO = [cr rcal - a sub ] + [a v 
+ [cr part + a s 2 uh ] , 



rtual 



sub] 



in which each bracket is separately finite.[] 

3. The results for the Tevatron are displayed in 
Figs. ||(a,b). For a Higgs mass between 100 and 150 GeV, 
the cross section varies between about 10 and 1 fb, the 
central value \x — ► /j,q = (2m t + Mh)/2 chosen for the 
renormalization and factorization scales. In NLO the 
theoretical prediction is remarkably stable with very lit- 
tle variation for /i between ~ fJ-o/3 and ~ 3/io, in contrast 
with the Born approximation for which the production 
cross section changes by more than a factor 2 within the 
same interval. Although at the Tevatron the cross section 
is strongly dominated by the qq annihilation channel for 
scales /i ~ /io, the proper study of the scale dependence 
requires inclusion of the gg, gq, and gq channels. If fi 
is chosen too low, large logarithmic corrections spoil the 
convergence of perturbation theory, and the NLO cross 
section would even turn negative for fi < /xo/5- 

As apparent from Fig. g(b), the K factor, K = 
cnlo/^lo with all quantities calculated consistently in 
lowest and next-to-leading order, varies from ~ 0.8 at the 
central scale fi — Ho to ~ 1.0 at the threshold scale fj, — 
2fj,Q- The small K factor can be understood intuitively 
in the fragmentation picture proposed in Ref.jlO]]. The 
average CM energy (vl) for the subprocess qq — > ttH 
at the Tevatron is about 650 GeV, i.e. sufficiently above 
the threshold region, so that the EHA of Ref. pOl can be 
used at least at a qualitative level, as confirmed earlier 
for e+e- -> ttH [g. For M% < m\ < (V$) 2 , the prob- 
ability for the hadronic process is decomposed into the 
product of probabilities for tt production and subsequent 
fragmentation t — ► t + H . As a result, the relative QCD 
corrections take the form S = 8[pp — ► tt] + 5[t — ► t + H]. 

With 8[pp -^qq-> it] a s /2ir Q and S[t -* t + H] ~ 

— 4a s /7r for small energies of the Higgs boson, the sum 
of the two terms 5 ~ — 9a s /2TT is negative and the K 
factor is predicted below unity in this limit. Integrating 
over the entire Higgs spectrum, the numerical evaluation 
yields ATeha ~ 0.7, which is nicely compatible with the 
result K ~ 0.8 of the full 0(a s ) calculation. 

Near threshold, Vs > 2mt + Mh, the QCD correc- 
tions are enhanced by Coulombic gluon exchange be- 
tween the top and antitop-quark in the final state. 
This Sommcrfeld rescattering correction |18| increases 
inversely proportional to the maximum t/t velocity 

/3 t max ~ \J(V~s - M H ) 2 - Am 2 t /2m t in the tt CM frame: 

<5coul(s) = C ■ na s /2 ■ (l//3 t > - C ■ 8a s /3/3J nax - H the 



tt pair is generated in a color-singlet state, the quark 
and anti-quark attract each other, and with C± = +4/3 



(3) 



As an independent cross check the phase-space slicing 
method has been applied to the subchannel qq —* ttH, which 
dominates at the Tevatron. The results obtained by the slic- 
ing and the subtraction techniques are in mutual agreement. 



the correction is positive. This leads to a strong en- 
hancement of the e + e~ — > ttH annihilation cross section 
near threshold ||. By contrast, if the tt pair is gen- 
erated in a color-octet state, the force is repulsive, and 
with Cg = —1/6 the correction is negative and relatively 
small. This applies to the dominant channel at the Teva- 
tron, qq — > ttH, which is mediated by s-channel color- 
octet gluon exchange. As a consequence, the destructive 
Coulomb interference term amplifies the reduction of the 
cross section. 
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FIG. 2. (a) The cross section for pp — > ttH + X at the 
Tevatron in LO and NLO approximation, with the renormal- 
ization and factorization scales set to /i = rat + Mh/2; (b) 
variation of the cross section with the renormalization and fac- 
torization scales for a fixed Higgs-boson mass Mb = 120 GeV. 

4. The improvement of the prediction for the cross sec- 
tion at the LHC , Figs. ||(a,b), is similarly striking. How- 
ever, the gluon initial states give rise to increased gluon 
radiative corrections [which will be improved by resum- 
mation techniques in the future]. For the central renor- 
malization/factorization scale (Iq we obtain K ~ 1.2, 



increasing to ~ 1.4 at the threshold value /z = 2^ - 
These values are nearly independent of Mh in the rel- 
evant Higgs mass range. 

The K factor at the LHC can be estimated in the 
fragmentation picture pp| , since the average subenergy 
(VI) ~ 830 GeV is relatively high at the LHC. With 
the dominant gg production channel, the sum of S[pp — > 
gg — ► tt] ~ +lla s /7r and 6[t—*t + H] ~ — 4a s /7T for the 
QCD corrections come now with opposite signs, but the 
positive correction to the tt production in the gg channel 
more than compensates the negative correction to the 
fragmentation, leading finally to 5 ~ +7a s /n. Taking 
into account the different renormalization and factoriza- 
tion scales, the estimate for the K factor in Ref. []10| is 
recognized compatible at the qualitative level with the 
full NLO result. 
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FIG. 3. Analysis of pp — » ttH + X at the LHC; (a) produc- 
tion cross section, and (b) renormalization/factorization-scale 
dependence (parameters as specified in the previous figure). 

5. Summary. The strong scale dependence of the Born 
cross sections in the reactions pp/pp — ► ttH , which pro- 



vide important search channels for the SM Higgs boson, 
requires the improvement by NLO QCD corrections. In 
agreement with a qualitative fragmentation picture, the 
K factor at the Tevatron is slightly below unity, i.e. vary- 
ing between ~ 0.8 and ~ 1.0 for renormalization and fac- 
torization scales /j, = fj,Q and 2 fig, with 2/xo denoting the 
threshold CM energy of the parton subprocesses. Sim- 
ilarly, the K factor varies between ~ 1.2 and ~ 1.4 for 
the same scale at the LHC. Most important, in contrast 
to the Born approximation, the predictions for the cross 
sections including NLO QCD corrections are stable when 
the renormalization and factorization scales are varied so 
that the improved cross sections can serve as a solid base 
for experimental analyses at the Tevatron and the LHC. 
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